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Abstract 
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show an interesting intermingling of the stripping reaction with elastic scattering of the reaction 
residue on the target. We also obtain doubly-differential distributions of the cross section on the 
parallel- and transverse- momentum variables. The distributions depend strongly on the value of the 
magnetic quantum number m. They will be of importance for evaluating acceptance corrections in 
experiments, and they lead to alignment with the possibility of anisotropic emission of subsequent 
gamma rays, an interesting spectroscopic tool. Experimental data for proton stripping of 8 B agree 
with our calculations. 
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I. INTRODUCTION 



Single-nucleon knockout reactions with heavy ions, at intermediate energies and in in- 
verse kinematics, have become a specific and quantitative tool for studying single-particle 
occupancies and correlation effects in the nuclear shell model, see the recent review The 
high sensitivity of the method has allowed measurements on rare radioactive species avail- 
able in intensities of less than one atom per second for the incident beam. The experiments 
observe reactions in which fast, mass A, projectiles with laboratory momentum = Pa/^ 
collide peripherally with a light nuclear target, typically 9 Be, producing residues with mass 
(A — l), in the following referred to as the core (c) of the assumed two-body system of core 
plus nucleon. In the laboratory system the momentum transferred to the core is 

A — l 

K,lab = ~\ kA - k A -i- (1) 



.4 

The final state of the target and that of the struck nucleon are not observed, but instead the 
energy of the final state of the residue can be identified by measuring coincidences with decay 
gamma-rays emitted in flight. Referred to the center-of-mass system of the projectile, the 
transferred momentum is k c . In the sudden approximation and for the stripping reaction, 
defined below, this must equal the momentum of the struck nucleon before the collision. 
The measured partial cross sections to individual final levels provide spectroscopic factors 
for the individual angular-momentum components j. In complete analogy to the use of 
angular distributions in transfer reactions, the orbital angular momentum I is in the knockout 
reactions revealed by the distributions of the quantity k c . These distributions are the subject 
of the present paper. 

The early interest in momentum distributions came from studies of nuclear halo states, 
for which the narrow momentum distributions in a qualitative way revealed the large spa- 
tial extension of the halo wave function. It was pointed out by Bertulani and McVoy 
that the longitudinal component of the momentum (taken along the beam or z direction) 
gave the most accurate information on the intrinsic properties of the halo and that it was 
insensitive to details of the collision and the size of the target. In contrast to this, the 
transverse distributions of the core are significantly broadened by diffractive effects and by 
Coulomb scattering. For experiments that observe the nucleon produced in elastic breakup, 
the transverse momentum is entirely dominated by diffractive effects, as illustrated by 
the angular distribution of the neutrons from the reaction 9 Be( 11 Be, 10 Be+n)X. In this case, 
the width of the transverse momentum distribution reflects essentially the size of the target. 
Experiments and theory for reactions of neutron halos have bee reviewed in ref. j4]. It was 
found that to understand the measured longitudinal momentum distributions it is necessary 
to take into account that a heavy-ion knockout reaction, being surface-dominated, can only 
sample the external part of the nucleon wave function. The magnitude of the reaction cross 
section is determined by the part of the wave function that is accessed, and the shape of the 
momentum distribution reflects the momentum content in this part. Calculations 0, 111 E| 
based on a sharp-surface strong-absorption ("black-disk") model could account for the ob- 
served longitudinal momentum distributions and also, approximately, for the absolute cross 
sections. This approach is confirmed in the, more accurate, present work, which extends the 
theory to include the general dependence of the differential cross section on the momentum 
vector. 

It is essential to note that the cross section for the production of a given final state 
of the residue has two contributions. The most important of the two, commonly referred 
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to as stripping or inelastic breakup, represents all events in which the removed nucleon 
reacts with and excites the target from its ground state. The second component, called 
diffractive or elastic breakup, represents the dissociation of the nucleon from the residue 
through their two-body interactions with the target, each being at most elastically scattered. 
These events result in the ejected nucleon being present in the forward beam with essentially 
beam velocity, and the target remaining in its ground state. These processes lead to different 
final states, they are incoherent, and their cross sections must be added in measurements 
where only the residue is observed. General expressions for the total and differential cross 
sections for the two components have been given by Hencken, Bertsch and Esbensen 0. 

In a subsequent development, the knockout method was extended to non-halo states 
pi. IqL fl~oL 11 . lil . 13l . For these, involving more deeply-bound nucleons, the one-nucleon 
stripping cross sections are much smaller than the free-nucleon reaction cross section on 
the same target; a ratio that gives a measure of how much the nucleon wave function is 
"shielded" from the target by the bulk of the core. This required a more elaborate theoretical 
treatment based on the elastic S- matrices S c and S n |l5 Jl6|] of the core and nucleon. For a 
general review of applications of this technique see ref. l|, which shows that very accurate 
theoretical single-particle cross sections are obtained in this way. However, the longitudinal 
momentum distributions have continued being calculated in the black-disk approximation. 
This simplification has been permissible because the assumed sharp surface of the target 
generates only transverse momentum components, which are integrated out in the final 
result. 

In the present paper we treat the three-dimensional momentum distribution of the core 
in stripping reactions at the same level of approximation as the single-particle cross sec- 
tions in Tostevin's calculations 0, Using S'-matrices from this work, we obtain identical 
single-particle cross sections after integration over momentum and summation over the m; 
substates. The results (i) test the reliability of the previously used sharp-surface approxima- 
tion, and (ii) demonstrate that the transverse momentum distributions are quantitatively 
and even qualitatively different from the parallel momentum distributions, and (iii) can serve 
to extract angular-momentum information from the angular distributions of the residues. 
Finally, (iv) the results are of importance for calculating acceptance corrections in experi- 
ments, and we evaluate the correlations between longitudinal and transverse components. 
After a presentation of the main features of the theoretical calculations, we give examples of 
the distributions obtained for different values of the nucleon separation energy S n and the 
orbital angular momentum /. Some essential numerical details are presented in an appendix. 



II. CALCULATION OF THE DIFFERENTIAL STRIPPING CROSS SECTION 

In this section we summarize the equations used in the calculation of stripping cross 
sections. The numerical calculations were done by using a Gaussian expansion method which 
is described in the appendix. We also show the equations used for the elastic scattering cross 
section, which will be important to interpret the numerical results obtained for the stripping 
cross sections, discussed in section 3. 

Following ref. 0], the cross section for the stripping reaction (c + n) + A — > c + X, 
where c corresponds to a specified final state of the core, is given by 



da. 



str 



d?k c (2tt) j 21 + 



d 2 K. 



\S n {K)\' 



d 3 r e" 



ik c .r 



S c (b c ) ipim (r) 



(2) 



3 




FIG. 1: Coordinates used in text. 



and where r = (p, z, . 



R r 



R c , so that 



p - b n | = Jp 2 + b 2 n - 2p b n cos i 



r 2 sin 9 + b 2 l — 2r sin # 6 n cos (<f> — </> n ) 



(3) 



with the symbols b c and b„ denoting two-dimensional vectors, which are the respective 
transverse components of R c and R n (see figure HJ). S c (S n ) are the S-matrices for the 
core+target and the neutron (or proton) -(-target scattering. 

The single-particle bound state wave functions for the subsystem (c + n), i.e. ipi m (r), in 
eq. (J2J) are specified by ipi m (r) = R[ (r) F; m (f), where Ri (r) is the radial wave function. It 
is not necessary to specify the total single-particle angular momentum j, since the assumed 
interaction is spin- independent, and the depth of the single-particle potential well is adjusted 
to reproduce the effective nucleon binding energy. 

The cross sections for the longitudinal momentum distributions are obtained by integrat- 
ing eq. (J2J) over the transverse component of k c , i.e. over k^, and using 



<i 2 kf exp — ik c .(p — p) = (2tt) 2 5 (p — p 



(4) 



One gets 



da. 



str 



1 1 



dk, 



1 - Jo 



1 Sn 



d 2 P \S c (b n )t 



(5) 



2tt2/ 

/oo 
dz exp [-ik z z\ ipi m (r) 
-oo 

where k z represents the longitudinal component of k c . 

For the transverse momentum distribution in cylindrical coordinates k± = ^Jk 2 , + k 2 , one 
uses in eq. (J2J) 

/oo 
dk z exp [— ik z (z — z')\ = 2ir5 (z — z') , (6) 
-oo 
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and the result is 



da K 



d 2 k l 



1 1 
2^2/ + 1 

TOO 



J2 J dz J d 2 p exp(-ikj.p) S c (b n )if), 



m, p 



(7) 



Sometimes it is convenient to describe the transverse momentum distributions in terms 
of the projection onto one of the Cartesian components of the transverse momentum. This 
can be obtained directly from eq. (|7J), i.e. 



da. 



str 



dkn 



dkr, 



da, 



str 



(8) 



The total stripping cross section can be obtained by integrating either eq. © or eq. (JJJ). 
For example, from eq. (JSJ), using eq. one obtains 



2tt 



C"str 



21 + 1 7 

d 3 r 



1-|^ (fen)! 1 



S'c ( Jr 2 sin 6 1 + b\ — 2r sin 6 1 b n cos i 



/ m (r) 



(9) 



which is the same as eq. (12) of ref. 0. 

In the Appendix A we show how the integrals in eqs. (0) and (0) can be evaluated 
numerically with use of Gaussian expansions of the core/nucleon+target S-matrices. 

The S-matrices have been obtained using the eikonal approximation for the wave func- 
tions. In this approximation the outgoing wave of a fragment, with wave number k, is given 
by 



(+) 



exp lik ■ r + 



hv 



dz'U, 



opt\ 



(10) 



The overlap of the incoming and outgoing wave function becomes 

*W|*W) = 5(6)fflq>(<q.r) , 



(11) 



where q = k' — k is the momentum transfer and S (b) is the scattering matrix given by 

1 



S (b) = exp [ix(b)} , 



with 



X(b) = ~ 



hv 



dz U opt (r 



(12) 



and Uopt(r) is the appropriate optical potential for the core+target and the neutron (or 
proton) +target scattering. In equation ()12j) \{b) is the eikonal phase, and r = \/b 2 + z 2 , 
where b is often interpreted as the impact parameter. This interpretation arises from a 
comparison of the results obtained with eikonal wavefunctions with those obtained with 
classical particles colliding at a given impact parameter b (l7j . Nonetheless, the eikonal 
wavefunction is a quantum scattering state and b is the transverse coordinate associated to 
it. Thus wave-mechanical effects, like smearing and interference, are accounted for properly. 

In the optical limit of the Glauber theory, the eikonal phase is obtained from the nuclear 
ground state densities and the nucleon-nucleon cross sections by the relation [I?} 



X(b) 



dq q P P (q) pt (q) Inn (q) Jo (qb) 



(13) 
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where Jo is the Bessel function of order zero, p Pt t (q) is the Fourier transform of the nu- 
clear densities of the projectile and target, and f^N (<?) is the high-energy nucleon-nucleon 
scattering amplitude at forward angles, which can be parametrized by |l8{ ] 

Inn (q) = ~f^~ aNN + aNN ) ex P (-(^NNq 2 ) • (14) 

In this equation <Jnn, &nn, and Pnn are parameters which fit the high-energy nucleon- 
nucleon scattering at forward angles. In eq. (|T3|) the quantities p p (q) and p t (q) are calcu- 
lated from the radial density distributions, usually [l| taken to be of Gaussian shapes for 
light nuclei, and of Fermi shapes for heavier nuclei with parameters taken from experiment. 
For cases where more accuracy is needed, it is possible to take the density distributions from 
Hartree-Fock calculations, as has been done in some recent work. This demonstrates that 
the theoretical uncertainty on the integral single-particle cross section (i.e. for a spectro- 
scopic factor of unity) is of the order of 5% for a halo state [19] and 15% for a very deeply 
bound I = 2 state (S n = 22 MeV) The precise choice of input parameters influences 
the absolute spectrosco pic factors, which have been the subject of many previous papers 
P. fsL Hi. 12l 3, Eil EH [H| to which we refer for numerical details. However, it means little 



for the shapes of the momentum distributions which are the focus of the present work, and 
it will not be discussed further here. 

For the Coulomb part of the optical potential the integral in eq. (I12J) diverges. One solves 
this by using \ — Xn + Xc, where \n is given by eq. (JT2J) without the Coulomb potential 
and writing the Coulomb eikonal phase, \c as 

X c(b) = 27 1 \n(kb) , (15) 

where 77 = Z p Z t e 2 /hv, Z p and Z t are the charges of projectile and target, respectively, v is 
their relative velocity, k their wavenumber in the center of mass system. Eq. (|15|) reproduces 
the exact Coulomb scattering am plit ude when used in the calculation of the elastic scattering 
with the eikonal approximation |17| : 

fc{6) = 2pv^f(6/2) 6XP H 7 ? ln [ sin W 2 )] + + 2*0 O } (16) 

where 0o — argT(l + irj/2). This is convenient for the numerical calculations since, as shown 
below, the elastic scattering amplitude can be written with the separate contribution of the 
Coulomb scattering amplitude included. Then, the remaining integral (the second term on 
the r.h.s. of eq. (|2Tj) below) converges rapidly for the scattering at forward angles. 

Although the Coulomb phase in eq. (fT5|) diverges at b = 0, this does not pose a real 
problem, since the strong absorption suppresses the scattering at small impact parameters. 
It is also easy to correct this expression to account for the finite charge distribution of the 
nucleus. For example, assuming a uniform charge distribution with radius R the Coulomb 
phase becomes 

Xc (b) = 2 V {e{b - R) \n(kb) + Q{R - b) [ln(kR) + ln(l + ^1 - b 2 /R 2 ) 

-^l-b 2 /R 2 ~\{l-b 2 /R 2 fl 2 ]} , (17) 

where B is the step function. This expression is finite for 6 = 0, contrary to eq. ()15p . If one 
assumes a Gaussian distribution of charge with radius R, appropriate for light nuclei, the 
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Parallel Momentum P z [MeV/c] 

FIG. 2: Longitudinal momentum distribution for the residue in the 9 Be( 11 Be, 10 Be 9S )X reaction 
at 60 MeV/nucleon as studied in the experiment of Aumann et al. The dashed curve is the 

cross section calculated in the transparent limit, re-scaled by a factor of 0.535. It has a high-energy 
component arising from the inner lobe of the wave function. The full drawn curve is the result of 
the present work. The black disk calculation (not shown) is essentially indistinguishable from this, 
except for the high-momentum tail. The dot-dashed curve shows that the weak shoulder is not 
present with an assumed nodeless 0s wave function with the same binding energy (0.504 keV). 



Coulomb phase becomes 

1 



Xcip) = 2r] 
where the error function E\ is denned as 



ln(M>) + |£i(6 2 AR 2 ) 



(18) 



E x (x) = I — dt . (19) 

J X t 

This phase also converges as b — > 0. In eq. (fT7j) R = R p + R t , while in eq. (jTSjl R = 
J Rp + Rt, where R p and R t are the respective projectile and target radius. The cost of 
using the expressions (fP7|) and (fT%J) is that the Coulomb scattering amplitude becomes more 
complicated than f)lfij) . Moreover, we have numerically attested that the elastic and inelastic 
scattering cross sections change very little by using eqs. (fT7j) or instead of eq. ()15|). 

In calculations involving stripping, the final state Coulomb interaction between the core 
and the target is taken into account by using the eikonal-Coulomb phase shift of ()15|) in the 
calculation of S c . However in the calculation of diffraction dissociation both S c and S n are 
calculated using the eikonal-Coulomb phase shift of (fT3j). 

The calculation of elastic scattering amplitudes using eikonal wave functions, eq. (fTUj). is 
very simple. They are given by |l7| 



POO 

U9)=ik dbb J (qb) {1-exp Mb)}} , (20) 
Jo 
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where q = 2/csin(#/2), and 9 is the scattering angle. The elastic scattering cross section is 
da e i/dVt = \ f e i(8)\ 2 . For numerical purposes, it is convenient to make use of the analytical 
formula, eq. (fTE|). for the Coulomb scattering amplitude. Thus, if one adds and subtracts 
the Coulomb amplitude, fc(0) in eq. (|2U|). one gets 

/•CO 

fM=fc{0)+ik / dbbJ (qb)exp[ixc(b)]{l-exp[i X N(b)]}- (21) 
Jo 

The advantage of using this formula is that the term 1 — exp [iXN{b)} becomes zero for 
impact parameters larger than the sum of the nuclear radii (grazing impact parameter). 
Thus, the integral needs to be performed only within a small range. In this formula, xc is 
given by eq. (|T3|l and fc(@) is given by eq. (fTfij) . with 

^-^C + fV-^-arctan-^-) , (22) 

where C = 0.5772156... is Euler's constant. 

The elastic cross section can be expressed in terms of the transverse momentum by using 
the relationships dVt ~ d 2 kj_/k 2 , and k± ~ q = 2k sin {0/2), valid for high-energy collisions. 



III. EXAMPLES OF MOMENTUM DISTRIBUTIONS 

We explore the consequences of the expressions developed in sect.|TTlby calculating mo- 
mentum distributions for selected cases. The longitudinal momentum distributions, corre- 
sponding to a projection on the beam {P z ) axis,turn out to be very close to those obtained 
in the simpler black-disk approximation 0, 0, . For this reason, a comparison with the 
numerous experimental data available is hardly necessary, but some references are given. 
A general discussion of longitudinal momentum distributions can be found in 0]. For the 
distribution projected to an axis perpendicular to the beam axis, the situation is different. 
We find in all cases an interesting intermingling of momentum components arising from 
stripping of the nucleon and from elastic scattering of the core fragment on the target. The 
latter mechanism is diffractive for light targets and Coulomb-dominated for heavy targets. 
There is very little useful experimental evidence on the transverse momentum distributions, 
as only the observation of coincident gamma rays can separate out the differential cross 
sections to individual final levels. (Even for halo nuclei, such as n Be and 15 C, 20-30% of 
the inclusive cross section goes to excited levels with different I values.) 

We have chosen to represent the momentum distributions graphically in the following 
way. Projections onto a single Cartesian coordinate are shown for one half axis only (the 
P z distribution is symmetric in the eikonal approximation). For cases where the magnetic 
quantum number m differs from zero, we have weighted the differential cross section with 
the multiplicity of 2, so that the sum over all m components gives the total cross section. 
The spatial momentum distribution does not depend on the azimuthal angle. It is therefore 
convenient to present it as a two-dimensional function of the parallel and the transverse 
momentum with the definition 

™±£±r, (23) 



dk±dk z d 2 k±dk z 

which normalizes to the total cross section when the integration is extended over the negative 
k 7 axis. 
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FIG. 3: Longitudinal momentum distribution for the residue in the 12 C( 8 B, 7 Be 9S )X reaction at 



76 MeV/nucleon as studied in the experiment of Enders et al. 14]. The proton binding energy 
for this proton halo state is 0.137 MeV. The dashed curves are the cross sections calculated in the 
transparent limit, re-scaled by a factor of 0.35. The narrower full drawn curves are the results 
of the present work. The thin curves show individual m = ±1 components with the sums given 
as thick lines. The black disk calculation (not shown) gives very similar results. Fig. 4 of ref. 
compares a similar (black-disk) calculation with experimental data. 



A. Longitudinal momentum distributions 



(24) 



The deuteron was the first halo system to be studied. Serber |21[ calculated the mo- 
mentum spectrum of neutrons from deuteron breakup on a light target in an approximation 
that amounts to replacing the core S-matrix S c (b c ) in eq. (J2J) by unity. This leads to the 
expression 

where <j v t is the proton-target reaction cross section, so that the stripping cross section (the 
term coined by Serber, see his comment in j23|) is that of a free proton, and the momentum 
distribution is the Fourier transform of the bound-state wave function. This is an excellent 
approximation for the deuteron, and it served later to explain qualitatively why early studies 
of stripping reactions of halo nuclei observed narrow transverse and longitudinal momentum 
distributions |2j,|24j. 

However, the transparent-core limit is never a good approximation for reactions of nuclei 
heavier than the deuteron. It overestimates the stripping cross section by a factor 2 for a 
pronounced halo system such as 11 Be, and by up to a factor of 25 for deeply bound nucleons. 
The momentum distributions are also modified in a significant way by the absorption given 
by the core S-matrix S c (b c ) as shown by the examples of the I = 11 Be neutron halo and 
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FIG. 4: Longitudinal momentum distributions for the reaction 9 Be( 34 Ar, 33 Ar(| + )) at 68.2 
MeV/nucleon. This I = 2 neutron-removal reaction leads to a final | + level bound by 18.42 
MeV. The solid curves represent the exact calculations and the dashed curves the sharp-cutoff 
approximation . 



the I =1 8 B proton halo shown in figs. [2] and H3 The underlying reason is that the knockout 
reaction of a heavy ion on an absorptive target, such as 9 Be, is surface dominated and 
samples only the momentum components there and in the tail of the wave function. The 
case of n Be in fig. [2] is good demonstration of this: The-high momentum components (the 
"shoulder" ) seen in the S c (b c ) = 1 approximation arise from the inner lobe in the wave 
function of the second s state and vanish almost completely in the knockout reaction, which 
cannot "see" the interior of the projectile. 

In the case of the 1=1 knockout on 8 B, shown in fig. El both m components are much 
narrower than in the transparent-core approximation, and the m=0 state is suppressed 
relative to the m=±l states. A simple geometrical argument |6] accounts for this: Projected 
onto the (x, y) plane, the two lobes of m=0 wave function are oriented along the z axis and 
shielded by the core from interacting (alone) with the target. 

The sharp-cutoff model 0, S has been a useful tool for discussing parallel-momentum 
distributions. In this approximation the S-matrices S c (b c ) and S n (b n ) (see fig.[TJ) are replaced 
by step functions with radii chosen to approximate the free core-target and nucleon-target 
cross sections by means of two parameters, a target radius and a minimum-impact parameter. 
In the following we have scaled the latter by a factor of typically 0.95 to get agreement with 
the, more accurate, total cross sections. This adjustment has essentially no influence on the 
calculated shapes. For the two halo cases discussed in figs.|21 and 01 there is essentially no 
difference from the more exact calculation. The same is true for the deeply bound s state 
shown later, in fig. [7| 

The case of a deeply bound state with 1=2 is shown in fig. The reaction 

9 Be( 34 Ar, 33 Ar(| + )) at 68.2 MeV/nucleon has a separation energy of 18.42 MeV and has 
been studied experimentally by Gade et al. . The solid curves are exact calculations and 
the dashed curves are obtained with the sharp-cutoff approximation. In this case there is a 
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FIG. 5: Transverse momentum distribution (full drawn) for the reaction 9 Be ( 11 Be, 10 Be gs ) at 
60 MeV/nucleon. The corresponding longitudinal momentum distribution (dashed) is shown for 
comparison. The dot-dashed curve is the calculated elastic cross section for the core. 

noticeable difference between the two approximations. 

Finally, we draw attention to an interesting effect that clearly is outside the scope of 
the present paper. The experiment on 34 Ar [25] shows an excess of intensity at the low- 
energy side of the 1=0 and 1=2 momentum distributions. This cannot be accounted for 
in the eikonal approximation as presented here, since it does not contain multistep time 
dependence, does not conserve energy, and yields symmetric line shapes. A similar low- 
momentum tail was observed in the I = momentum distributions from neutron knockout 
on the n Be and 15 C halo states [26|. In this case, it was interpreted as arising from the 
diffractive reaction channel and could be quantitatively accounted for in a coupled-channels 
calculation. However, diffraction dissociation is a small contribution to the cross section 
in the case of 34 Ar and is hardly the explanation here. We have no explanation for this 
asymmetry, which has been seen also in other cases, and which remains an interesting 
problem for future investigations. 



B. Transverse momentum distributions 

Transverse momentum distributions depend more strongly on the details of the nucleus- 
nucleus interaction than do the longitudinal momentum distributions 0. The nuclear size, 
the diffuseness of the nuclear matter distribution, and the core-target Coulomb repulsion 
all contribute to the transverse momentum distributions. At large impact parameters the 
Coulomb force still has a strong influence, especially for heavy targets, which for halo systems 
make Coulomb dissociation the dominant breakup channel. 

In figure we compare the longitudinal and transverse momentum distributions for the 
reaction 9 Be ( 11 Be, 10 Be) at 60 MeV/nucleon. The transverse momentum distributions (solid 
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FIG. 6: Transverse momentum distribution (full drawn) for the reaction 9 Be( 15 C, C gs ) at 103 
MeV/nucleon. The parallel-momentum distribution (dashed) has been studied experimentally 
T^ . 2f|. The influence of the Coulomb interaction is small, but visible with the log scale, for this 



target-core combination. The curves labelled "no Coulomb" and "(-Z)" represent a calculation 
with the target charge set to zero and to -4, respectively. 



curve) has for small momenta the same shape as the longitudinal one (dashed), but is lower 
in intensity by about 10%. The missing cross section shows up as a broad distribution 
with an oscillatory pattern, which we interpret as elastic scattering of the core fragment 
simultaneously with the stripping of the neutron. We demonstrate this by calculating the 
core-target elastic cross section, which shows a similar pattern with the minima characteristic 
of Fraunhofer scattering. The minima are not in the same place, however, because of the 
condition that the neutron be absorbed in the stripping reaction. A similar broad component 
appears in the transverse momentum distribution of the core fragment in coupled-channels 
calculations (ref. |26j] and J. A. Tostevin, personal communication). 

For the case of 11 Be, bound by 0.504 MeV, the large size of the halo allows most absorption 
processes to occur without simultaneous elastic scattering of the core. For more bound states 
the two mechanisms become increasingly intermingled. Already in the stripping of the halo 
neutron of 15 C, bound by 1.218 MeV, the tail from core scattering is much stronger, as 
seen in fig.® and for the 1=0 neutron knockout from 34 Ar, bound by 17.06 MeV, the two 
mechanisms are no longer distinguishable, as seen in fig-El The case of knockout of a 
deeply bound 1=2 neutron, also from 34 Ar, is shown in fig.|HJ which again shows a transverse 
momentum distribution that is broader than the longitudinal distribution shown in fig. HI 
Note the different momentum dependence of the individual m components, which can easily 
be understood from geometrical properties of the spherical harmonics 0. 

The longitudinal momentum distributions are not affected by elastic scattering of the 
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FIG. 7: Transverse momentum distribution (full drawn) for the 1=0 knockout reaction 
9 Be( 34 Ar, 33 Ar(^ + ) at 68 MeV/nucleon and with a neutron separation energy of 17.06 MeV j^jj. 
The corresponding longitudinal momentum distribution (dashed) is essentially indistinguishable 
from the calculation in the black-disk approximation (dot-dashed). 



core fragment. The basic reason for this is that the forces acting on the core fragment are 
along the z-direction and reverse sign at the origin of the z axis, thus leading to a null effect 
provided the angular deflection of the core fragment is small. For the same reason, the 
longitudinal momentum distributions are also relatively insensitive to Coulomb effects. 

We now examine how the Coulomb force affects the transverse momentum distributions. 
In the case of a neutron halo, it pushes the core away from the target along the transverse 
direction. The resulting distributions for the case of a light target are here illustrated by 
carbon (Z—6) on beryllium (Z=4), shown in fig.El The effect is small and would hardly be 
measurable. We note, however, that leaving out the Coulomb interaction gives a broader 
momentum distribution, suggestive of an interference effect. There is an interesting way to 
explore this. 

It is well known that the higher order Coulomb effects, e.g. to second-order in perturba- 
tion theory, carry information on the absolute sign of the charge of the particles. In atomic 
physics this is often referred to as the Andersen-Barkas effect, which reflects contributions 
of odd powers in the projectile charge on stopping powers and ionization probabilities. In 
nuclear physics this idea has been explored in the study of dynamical effects in Coulomb 
dissociation j22|. We have carried out calculations reversing the sign of the target's charge 
but keeping all nuclear interactions constant. Fig.El shows that a hypothetical core-target 
Coulomb attraction, as could be expected, gives an even broader distribution. The reason 
for this is that the elastic scattering at large transverse momenta is dominated by the far-side 
scattering contribution, in which the nuclei pass by close to each other. For a negatively- 
charged target both nuclear and Coulomb forces interfere constructively, leading to larger 
deflection angles. 

The Coulomb effect becomes important for heavier targets. To illustrate this, we compare 
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FIG. 8: Transverse momentum distribution (full drawn) for the 1=2 knockout reaction 
9 Be ( 34 Ar, 33 Ar) at 68 MeV/nucleon, corresponding to the | + excited level in the core with an ef- 
fective neutron separation energy of 18.43 MeV. The longitudinal momentum distribution (dashed) 
is shown for comparison. 

in figure El the transverse momentum distributions for the ( 11 Be, 10 Be) reaction on silicon 
(upper panel) and on lead (lower panel) targets. One sees that the momentum distributions 
are somewhat distorted by the Coulomb field for a silicon target, while for a lead target 
the effect is huge and the outcome is completely dominated by Coulomb repulsion. The 
transverse momentum transfer to the core fragments after the stripping can be roughly 
estimated by the relation AP C = 2Z c Zxe 2 /bv , where h is the impact parameter and v is 
the projectile velocity. For the reaction on a lead target, the minimum impact parameter 
is b ~ 10 fm, and with v ~ 0.4 c, we obtain AP C ~ 220 MeV/c. This agrees well with 
the upper limit of the calculated distribution; the broad peak (lower panel of figure EJ) is 
approximately one half of that value. 

Finally, we remind the reader that the present paper is an analysis of the stripping re- 
action, also called inelastic breakup. For experiments that only observe the core fragment, 
there is also an incoherent contribution from elastic breakup, which is given by the coher- 
ent sum of the contributions from diffraction dissociation and Coulomb dissociation. The 
latter becomes dominant for the reaction of 11 Be on lead and will lead to momentum 
distributions that, in principle, are different from those shown in fig. 01 

IV. CALCULATION OF THE DOUBLE-DIFFERENTIAL STRIPPING CROSS 
SECTION AND COMPARISONS WITH EXPERIMENT 

A. Alignment of the reaction residues 

Except for the special case in which the wave function factorizes in Cartesian coordinates, 
the longitudinal and transverse momentum components will be coupled. This is of primary 
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FIG. 9: Transverse momentum distributions of 10 Be fragments for the reaction ( 11 Be, 10 Be) at 
80 MeV/nucleon. The solid (dotted) curves include (do not include) the core-target Coulomb 
interaction. The upper panel (a) is for Si targets, while the lower panel (b) is for Pb targets. 
Changing the sign of the charge of the target yields the dashed curves. 

importance for the analysis of experiments in which the detection system limits the accep- 
tance of events to a certain volume in momentum space, a problem that until now has been 
treated in simple approximations. The selection in (P z , P±) space may also in special cases 
provide an additional spectroscopic tool. In order to illustrate the phenomenon, we plot in 
fig- El the double differential cross section d 2 a/dP z dP± for the reaction of a deeply bound 
1=2 state calculated from eq. ()23|) . The momentum distributions for the I = and the three 
I = 2 states are shown separately. For the case 1=2, the distributions for the three m states 
are peaked in different regions of the P z -P± plane. The maps of fig. El are clearly necessary 
input to an accurate calculation of the experimental acceptance. 

We have already pointed out that the reactions favor the formation of residues in states 
with the maximum absolute value of the magnetic quantum number m = ±Z (with the 
quantization axis taken to be the beam direction). This (possible) alignment effect implies a 
(possible) anisotropic emission (in the center-of-mass system) of gamma rays emitted after 
the reaction. This is a tool for identifying multipolarities of the gamma transitions and 
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FIG. 10: Contour plots for the 1=0,2 knockout reactions 9 Be ( 34 Ar, 33 Ar) X at 68.2 MeV/nucleon. 
The contour lines are equidistant. The projections on one Cartesian coordinate axis were illustrated 
in figs-EJEJ and|Hl The absolute values of d 2 a/dk z dk± (in fm 2 / (MeV/c) 2 ) are given by the numbers 
in the contour plots, which are to be multiplied by the factors shown in each panel. 



spin sequences in the product nucleus. The fact that for different values of m, the main 
contributing cross sections are located in different areas of the (P z , P±) map suggests that 
the alignment effect can be enhanced by making cuts in the momentum components. A 
theoretical example of such an application has been given in fig. 12 of ref. £0], which shows 
the calculated angular distributions for two different spin sequences. For the example given 
in fig. ^2 the ^=2 fraction in the reaction residues is ( "set 1" ) 58% as compared with 40% 
for a population with statistical weights. Selecting reactions with P z values between -50 and 
50 MeV/c ("set 2") increases the m=2 fraction to 85% and reduces the count rate to one 
half, corresponding to a net gain in sensitivity. Limiting in addition the values of P± to 
values between 85 and 165 MeV/c gives only a marginal improvement: The m=2 fraction 
increases to 88% but at the cost of reducing the intensity by another factor of one half. 

As an example of a possible application of the alignment effect, we show in fig. ITT1 angular 
distributions of gamma rays in the center-of-mass system calculated from the expressions 
given by Yamazaki [28|. The example is for an assumed | + to | + cascade as in 33 Ar and 
for various assumed multipoles and mixing ratios. The alignment parameters correspond to 
"set 2", which gives results that are already approaching pure m=±2 alignment. Without 
the cut on P z the anisotropy would be only half as big. 
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FIG. 11: Gamma-ray angular distributions for a | to j transition with alignment parameters 
corresponding to "set 2" (in the text) obtained by a central cut on the parallel momentum. Calcu- 
lations are shown for the multipolarities E2, Ml, and for mixed transitions with E2/M1 amplitudes 
of ±1.0. 



B. Comparisons with experiment 

The present work deals with an extension of a theory jj| 0] that has been tested in 
numerous experiments. As indicated in subsect. flll Al the longitudinal momentum distri- 
butions are well in hand, both experimentally and theoretically. There is little reason to 
doubt that the same is the case for the transverse-momentum distributions calculated here. 
However, there are few good data sets to compare with. Essentially all measured transverse 
distributions that we are aware of are superpositions of contributions from several / values 
(see the large set of inclusive momentum spectra reported by Sauvan et al. [29]). For these, 
coincidence measurements with gamma rays would be required to provide detailed test of 
our calculations. However, the proton knockout on 8 B is an exception to this. 

For 8 B, the inclusive transverse-momentum distribution in the knockout of the halo proton 
has been measured in the experiment by Kelley et al. 3(|. In spite of the absence of 



gamma-ray coincidence data, this is a favorable case because the ground-state cross section 
dominates and because the approximately 15% branch to the excited level also has /=1 and 
has an almost identical shape. Fig.1121 shows that our calculation is in excellent agreement 
with the data, which it reproduces over two orders of magnitude in cross section. There are 
several measurements of the longitudinal- momentum distribution for 8 B proton knockout; 
fig. 4 of ref. 0] presents a comparison of high-energy data with a black-disk calculation, 
essentially indistinguishable from fig.|3]of the present work. 

Our second comparison with experimental data illustrates how the calculations in the 
present work may be applied toward clarifying a much more complicated issue. The reactions 
and structure of the two-neutron halo nucleus n Li have attracted much interest. It is 
a Borromean system in the sense that although the three-body system, consisting of 9 Li 
and two neutrons, forms a bound state, none of the possible two-body subsystems have 
bound states. Hence the stability of 11 Li is brought about by the interplay of the core- 
neutron and the neutron-neutron interactions, which must lead to a strongly correlated 
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FIG. 12: Inclusive transverse-momentum distribution for the residue in the 9 Be( 8 B, 7 Be 3S )X re- 
action measured at 41 MeV/nucleon [3(J. The theoretical calculation (full drawn) is based on 
the same parameter set as fig-EH and the only adjustable parameter is the (arbitrary) scale. The 
binding energy for this proton halo state is 0.137 MeV. The dot-dashed and dashed curves are the 
individual contributions of the 7 Be ground state and first-excited state, respectively. The angular 
resolution in the experiment broadens the data by approximately 4%. This has not been included 
in the theoretical curves. 



wave function with the two neutrons spatially close together. Barranco et al. |^l( found 
that attempts to understand n Li via breakup reactions are made difficult by important final- 
state interactions, so that the primary mechanism is removal of a single neutron followed by 
the decay in flight of the unstable 10 Li. The slow neutron emerging from this decay carries 
no direct information on the n Li structure; its energy spectrum reflects properties of both 
the initial and final state. 

Simon et al. showed that this problem can be circumvented by reconstructing the 
momentum vector of the 10 Li intermediate from the observed momenta of 9 Li and a neutron 
following a reaction on a carbon target. In an experiment at 287 MeV/nucleon they obtained 
the best resolution for the projected transverse- momentum component and arrived at the 
spectrum shown in fig. 113) T he momentum resolution (full width at half maximum) was 
approximately 55 MeV/c|33| . Their analysis was based on a simplified model using analyt- 
ical approximations to the momentum distributions and found the spectrum to require an 
(lsi/ 2 ) 2 component of (45±10)%, the rest being (0pi/ 2 ) 2 . They obtained further confirma- 
tion of this interpretation by measuring the angular distribution of the neutron relative to 
a quantization axis taken along the recoil direction of the 10 Li composite. This distribution 
showed the forward-backward asymmetry characteristic of interference between final states 
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FIG. 13: Inclusive transverse- momentum distribution for the residue in the inclusive 
12 C( 11 Li, 10 Li)X reaction measured at 287 MeV/nucleon [32| . The theoretical result (full drawn 
curve) is the adjusted sum of three angular-momentum components as described in the text. 

of opposite parity. 

For the calculation here we approximate the 11 Li ground state by an inert 9 Li core coupled 
to a neutron pair in a mixture of (lsi/2) 2 , (Opi/2) 2 , and (Og^) 2 states. We assume a two- 
body model, thus neglecting interference effects, and adjust the single-particle wave functions 
to reproduce the effective neutron separation energies. The two-neutron separation energy 
is 0.3 MeV and to this comes the average excitation energies of the three states in 10 Li. 
From the systematics in fig. 6 of (34| augmented with data for the 0d 5 / 2 neutron states 
in neighboring nuclei we estimate 10 Li average excitation energies to 0.2, 0.5, and 1.5 MeV 
for the three single-particle states, respectively. In the spirit of the sudden approximation, 
we take these to be the center of gravity for the multiplets formed by the coupling to the 
9 Li spin and arrive at effective neutron-separation energies of 0.5, 0.8, and 1.8 MeV, still 
in the same order, which determine the relative core-particle wave functions of the system 
10 Li+n. The sudden approximation suggests that the same wave functions should be used for 
calculating the S matrices (profile functions) for 10 Li treating it as a neutron+ 9 Li composite 
as described in ref. |9(. The three spectroscopic factors are then obtained from an adjustment 
to the differential cross section of fig.fT3l 

The calculated transverse-momentum distributions for the assumed three components 
were folded with the experimental resolution. The sum of the resulting distributions, 
weighted with the unknown spectroscopic factors, were adjusted in a \ 2 analysis to the 
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experimental result as shown in fig.1131 Per degree of freedom, the goodness of the fit 
is x 2 l v — 1-38 for absolute spectroscopic factors of 0.98(6) (lsi/ 2 ), 1.91(16) (0pi/ 2 ), and 
0.12(10) (Oo?5/2). In view of the many approximations made in this analysis, it is probably 
satisfactory that the the sum of the spectroscopic factors is 3.0, where the sum-rule value 
is 2. Re-normalized to the sum, the relative contributions are 33(2)%, 64(5)%, and 4(3)% 
in approximate agreement with the analysis of Simon et al. [32] cited above, who obtained 
45% for the lsi/ 2 state. The errors cited do not include systematic contributions arising 
from the theory or contained in the experimental data. It is interesting to compare with the 
neighboring nucleus 12 Be ^H, also with 8 neutrons, for which the 0<i5/2 contribution to the 
one-neutron knockout reaction appears to be much larger, of the order of 50%. 

V. CONCLUDING REMARKS 

The present work has extended the theory of one-nucleon stripping reactions (inelastic 
breakup) to cover two-dimensional momentum distributions of the reaction residue with the 
use of realistic profile functions for the core-target and nucleon-target interactions. The 
parallel-momentum distributions, projections on the beam direction P z , are not appreciably 
different from those obtained in the "black disk" approximation used in earlier work. On 
the other hand, the projections on the transverse direction, here referred to as P y , are very 
different from the projections on P z . For halo systems they show a weak component, which 
represents an additional mechanism in which the residue has scattered elastically on the 
target. For more strongly bound systems the two processes, stripping and elastic scattering, 
become inseparable, and the distributions on P y are broader than those on P z . 

The doubly-differential distribution of the cross section on the parallel- and transverse- 
momentum variables (P Z ,P±) has for a given angular momentum / > 1 a characteristic 
behavior for the different components of the magnetic quantum number m. This will be of 
importance for evaluating acceptance corrections in experiments. The m dependence of the 
cross sections leads to alignment of the residues and consequently to anisotropic emission 
of subsequent gamma rays. This effect can be exploited for identifying spin sequences 
and gamma-ray multipolarities as illustrated in the example in fig. ^2 A measurement of 
the transverse momentum distribution in stripping of 8 B on a light target is in excellent 
agreement with our calculation, which reproduces the experimental result accurately over 
two orders of magnitude in cross section. We have also attempted a more tenuous application 
of the theory to the complex case of neutron knockout on the two-neutron halo nucleus 11 Li. 
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VI. APPENDIX: GAUSSIAN EXPANSION METHOD FOR STRIPPING REAC- 
TIONS 



Using the explicit form of the spherical harmonics 
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Part of the integral in (J2J) is 
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To simplify the calculations we can express S c (b c ) as an expansion in terms of integrable 
functions. The S-matrices can be well described by the expansion 
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S c (b c ) = J2 a i ex P -bl/fi 
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(28) 



Good fits for realistic S-matrices were obtained with N = 20, i.e. with 20 complex coefficients 
otj and Rl = 10 — 20 fm, depending on the size of the system. Since the real part of the 
S-matrices has the usual behavior of S c (b c ) ~ for b c R, and S c (b c ) ~ 1 for b c ^> R, 
where R is a generic nuclear size, one of the coefficients of the expansion in eq. (J28|) is 
(Xj = 1, and Bj = oo, which we take as the j = term in the expansion. 
We now use the sum (|25j) and the equation 
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(29) 



valid for any complex z, in equation ()27jl. 

The integration in eq. (|27|) will then involve functions of the form 
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where r = (p,z), cos 8 = z/r = z/y/p 2 + z' 2 . Then 

Fim, j{k±,k z ,b n ) = Ci m aj exp -fi£//3? J dp p exp -p 2 /p 2 
x J dz exp [— ik z z) Ri{r) Pi m (cos9) 

e m ^ exp {— ikj_p cos (0^ — 0)} exp 

Using the expansion (J2HJ), we can write 
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Upon squaring eq. inserting in eq. (j2J), and integrating over 0^ and n , we can use 

J d(j) e^e"^ = 2tt(5 Pi p/ and we get 
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Since 



i I \ ( 00 ) % J ry ( liOC } 



where I a (x) is the modified Bessel function, one gets 
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where an irrelevant phase i m p was dropped off, as only the absolute value of Ai mp enters 
eq. flU}. 

The first term of the equation (jHHJ), with (3j = oo and ctj = 1 can be calculated using 

la (0) = 5 a . 

Using the integral 
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in eq. ()36|) one gets for the longitudinal momentum distribution 
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Using the integral of eq. © in eq. (|36J) one gets for the transverse momentum distribution 
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The formulas above are also useful to check the quality of the Gaussian fit, eq. (J28j) . to 
obtain the momentum distributions. One can compare the direct numerical integrations 
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using eq. (f4"Uj) with 
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